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Abstract 

In mammals, the circadian clock allows them to anticipate and adapt physiology around the 24 hours. Conversely, 
metabolism and food consumption regulate the internal clocl<, pointing the existence of an intricate relationship between 
nutrient state and circadian homeostasis that is far from being understood. The Sterol Regulatory Element Binding Protein 1 
(SREBPl) is a key regulator of lipid homeostasis. Hepatic SREBPl function is influenced by the nutrient-response cycle, but 
also by the circadian machinery. To systematically understand how the interplay of circadian clock and nutrient-driven 
rhythm regulates SREBPl activity, we evaluated the genome-wide binding of SREBPl to its targets throughout the day in 
C57BL/6 mice. The recruitment of SREBPl to the DNA showed a highly circadian behaviour, with a maximum during the fed 
status. However, the temporal expression of SREBPl targets was not always synchronized with its binding pattern. In 
particular, different expression phases were observed for SREBPl target genes depending on their function, suggesting the 
involvement of other transcription factors in their regulation. Binding sites for Hepatocyte Nuclear Factor 4 (HNF4) were 
specifically enriched in the close proximity of SREBPl peaks of genes, whose expression was shifted by about 8 hours with 
respect to SREBPl binding. Thus, the cross-talk between hepatic HNF4 and SREBPl may underlie the expression timing of 
this subgroup of SREBPl targets. Interestingly, the proper temporal expression profile of these genes was dramatically 
changed in Bman~'~ mice upon time-restricted feeding, for which a rhythmic, but slightly delayed, binding of SREBPl was 
maintained. Collectively, our results show that besides the nutrient-driven regulation of SREBPl nuclear translocation, a 
second layer of modulation of SREBPl transcriptional activity, strongly dependent from the circadian clock, exists. This 
system allows us to fine tune the expression timing of SREBPl target genes, thus helping to temporally separate the 
different physiological processes in which these genes are involved. 
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Introduction 

Mammals possess an internal circadian clock which allows them 
to anticipate and adapt to daily environmental changes [1]. The 
molecular mechanism underlying the cell-autonomous circadian 
rhythms relies on a network of feedback loops in which BMAL 1 , 
CLOCK, Neuronal PAS domain (NPAS) protein 2 and Retinoic 
acid receptor-related Orphan Receptor (ROR) proteins act as 
transcriptional activators and period homolog proteins (PERI, 2 
and 3), cryptochromes (CRYl and 2) and REV-ERBs function as 
inhibitors producing the self-sustained oscillating production of 
tiieir target genes, including themselves. At central level, the 
expression of clock genes is dictated by a pacemaker localized in 
the hypothalamic suprachiasmatic nucleus (SCN) which synchro- 
nizes the phase in nearly all body cells. However, in peripheral 



organs such as the liver, oscillations are also entrained by the 
feeding and fasting cycle [2,3,4]. This sophisticated regulatory 
system contributes to coordinate many physiological processes, 
such as sleep-wake cycles, locomotor activity, body temperature, 
hormone secretion and energy metabolism that all display 
circadian rhythms. In particular, the importance of the connection 
between circadian clock and metabolism regulation is emerging. 
Epidemiological studies have shown an increased incidence of 
obesity, diabetes, and cardiovascular disease, in addition to certain 
cancers and inflammatory disorders in night workers [5—7]. 
Accordingly, in genetic mouse models the disruption of the clock 
alters metabolic homeostasis at different levels (reviewed in [8]), 
suggesting a still unresolved relationship between nutrient state 
and circadian homeostasis. 
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Author Summary 

CIrcadian rhythmicity is part of our innate behavior and 
controls many physiological processes, such as sleeping 
and waking, activity, neurotransmitter production and a 
number of metabolic pathways. In mammals, the central 
circadian pacemaker in the hypothalamus is entrained on a 
daily basis by environmental cues (i.e. light), thus setting 
the period length and synchronizing the rhythms of all 
cells in the body. In the last decades, numerous investi- 
gations have highlighted the importance of the internal 
timekeeping mechanism for maintenance of organism 
health and longevity. Indeed, the reciprocal regulation of 
circadian clock and metabolism is now commonly accept- 
ed, although still poorly understood at the molecular level. 
Our global analysis of DNA binding along the day of Sterol 
Regulatory Element Binding Protein 1 (SREBPl), a key 
regulator of lipid biosynthesis, represents the first tool to 
comprehensively explore how its activity is connected to 
circadian-driven regulatory events. We show that the 
regulation of SREBPl action by nutrients relies mainly on 
the control of its subcellular localization, while the 
circadian clock influences the promoter specific activity 
of SREBPl within the nucleus. Furthermore, we identify the 
Hepatocyte Nuclear Factor 4 (HNF4) as a putative player in 
the cross-talk between molecular clock and metabolic 
regulation. 

The Sterol Regulatory Element Binding Protein 1 (SREBPl), a 
basic Helix-Loop-Helix-Leucine Zipper (bHLH-LZ) transcription 
factor, plays a key role in the regulation of lipid biosynthesis, which 
is one of the most feeding-related function in the liver [9]. 
SREBPl is synthesized as an inactive precursor, anchored to the 
ER-membrane and its N-terminal fragment is released into the 
nucleus after proteolytic cleavage in response to cholesterol 
depletion [10] or to activation of the insulin signalling pathway 
[11,12]. Two SREBPl isoforms, la and Ic, are obtained through 
alternative splicing of the same gene [13,14]. The liver expresses 
mostly the SREBPl c isoform that mediates the insulin-driven 
lipogenic activity [15]. Besides being under the control of the 
feeding-fasting cycle, SREBPl translocation to the nucleus is also 
influenced by one of the master clock regulators, REV-ERBa [16]. 
Nevertheless, in absence of a functional clock, such as in ctyl ' ; 
cry2~^~ mice, a normal expression pattern of several SREBPl 
target genes can be restored by an imposed rhythmic food intake 
[4] , suggesting a dominant role of the feeding-fasting cycle in the 
regulation of SREBPl. 

To systematically understand how the interplay of circadian 
clock and nutrient-driven rhythm regulate SREBPl activity, we 
evaluated the genome-wide binding of SREBP 1 to its targets along 
the day in wild-type mice. Our results define SREBPl binding 
pattern in the physiological context of both rhythmic food 
absorption and circadian rhythm and they give the first tools to 
comprehensively explore how SREBPl activity is connected to 
circadian-driven regulatory events. 

Results 

SREBPl binding to DNA is rhythmic 

To evaluate the genome- wide dynamics of SREBP 1 binding to 
its target sites in a physiological context, we prepared liver 
chromatin from C57BL/6 mice, collecting samples each 4 hours 
during one day (see Material and Methods). ChlP-seq with an 
antibody that recognizes both SREBPl isoforms was performed at 
each time point. The SREBPl antibody was tested extensively 



(Figure SI) and has also been used in previous studies [17]. We 
obtained an average of 38 millions sequence reads by time point 
by ultra-high-throughput sequencing (Table SI). The mapping 
allowed the identification of 448 bona fide SREBP 1 binding peaks, 
above the background. As shown in Figure lA, the binding of 
SREBP 1 is overall oscillatory, with the maximum for most of the 
sites at Zeitgeber Times (ZT) 1 4 or 18 (light is on at ZTO and is off 
at ZT12). To systematically evaluate the rhythmicity of SREBPl 
recruitment to its targets, a cosine function was fitted to the 
temporal profile of the binding (see Material and Methods). This 
allowed to calculate, for each peak, the binding phase and the 
amplitude of the oscillation together with its associated P-value, as 
exemplified for the two sites found on the Srebpl gene itself 
(Figure IB). 53% of SREBPl binding sites were found to be 
rhythmic (P<0.1 for the amplitude). Four clusters of targets were 
clearly distinguishable based on binding kinetics (Figure lA), the 
first with a phase distributed around ZT15-ZT17, whereas the 
other ones with the phase peaking around ZT11-ZT12, as 
determined by cosine function (Figure IC). The observed kinetics 
of SREBP 1 binding, especially for cluster A peaks, was consistent 
with its gene expression and nuclear localization, (Figure ID, IE 
and IF). Collectively, these results show that the activity of 
SREBPl oscillates with a pronounced circadian rhythm, in 
agreement with the previously reported daily variations of its 
RNA and protein levels [16,18-20]. 

Characterization of SREBPl binding sites 

SREBPl binding sites identified in this study are grouped in 
four clusters with a slightiy shifted phase. To better investigate the 
features of these sites we calculated the number of nucleotides 
spanned by each peak and found that sites belonging to cluster A 
(236 out of 448) were narrow, with a typical length of about 200 
nucleotides (Figure 2A). These peaks were also closer to the nearest 
annotated transcription start site (TSS) than peaks belonging to 
clusters B, C or D, whose distance to the nearest annotated TSS 
roughly matches randomly picked genomic locations (Figure 2B). 
Moreover, the amplitude of the binding oscillation along the 
24 hours was greater for cluster A peaks (Figure 2C). These 
observations suggest that cluster A sites may be more relevant in 
the regulation of transcription mediated by SREBPl. In agree- 
ment with this hypothesis, a MEME [21] motif search analysis 
clearly identified the canonical SREBP 1 consensus motif in more 
than 60% of the sites belonging to cluster A (Figure 2D), but only 
in 6% of the sites assigned to clusters B, C and D. Within cluster A, 
the MEME analysis also showed that motifs for SPl and NFY, two 
transcription factors (TFs) known to cooperate with SREBPl to 
regulate the transcription of its target genes, were overrepresented 
[17]. The consensus motif for the Hepatic Nuclear Factor 4 
(HNF4) was also identified in 6 1 out of the 236 cluster A sites. The 
discovered motifs were enriched in cluster A peaks with an 
empirical p-value<0.001, as shown in Table S2. In addition, in 
the regions belonging to cluster B, C and D, we determined only 
several highly repetitive sequences as top-scoring motifs (for 
example AC AC AC AC A in 73 sites out of 212) that could not be 
associated to any known consensus motif for TFs. This result 
suggests that SPl, NFY and HNF4 may participate to SREBP 1- 
mediated transcriptional regulation and further supports the 
functional importance of SREBPl binding sites assigned to cluster 
A. Thus, we opted to focus the following analyses on these regions, 
although we cannot exclude that the other sites might contribute 
to mediate SREBPl activity in mouse liver potentially through 
genome loops. 

To explore the cellular processes that are regulated by SREBP 1 
along the day, we annotated each site with the nearest Ensembl 
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Figure 1. Dynamics of SREBP1 binding. (A) C57BL/6 mice were fed only during the night (ZT12-ZT24) for one week before collecting liver every 
4 hours for one day. Chromatin from 5 mice was pooled at each time point and ChIP with an antibody against SREBP1 was performed. Peaks were 
positioned where the signal for SREBP1 was at least a four-fold in comparison to the input signal in at least one time point. The heat-map represents 
SREBP1 binding to all its targets along the time. Hierarchical clustering was done using Pearson correlation scores and identified four major clusters 
(A, B, C and D). The color scale is indicated below. In the column on the right, black lines indicate that Pol II was detected in the same site as SREBPl in 
at least one time point, as assessed in our previous ChlP-seq data set [24], (B) Two SREBPl binding sites were identified in the proximity of Srebpl 
gene, at a distance of -30 and 4-299 nucleotides from the Srebplc and Srebpla TSS, respectively. Graphs represent the fitting to a cosine function of 
experimental data obtained on these peaks (black dots), in order to calculate the phase of the binding (dashed line), its interval of confidence (dotted 
lines) and the associated P-value. (C) Histogram of binding phase frequency in clusters A, B and C, for peaks with a P-value of the amplitude <0.1. 
None of the peaks belonging to cluster D met this requirement. (D) mRNA expression of Srebplc was evaluated by qPCR in livers from C57BL/6 mice 
at the indicated ZT time (n = 5). Data are normalized using 36B4 as housekeeping gene. (E) Hepatic nuclear extracts from C57BL/6 mice were 
subjected to western blot analysis to detect the nuclear SREBPl . U2AF was used as loading control. Each sample is a pool of 5 livers. (F) Quantification 
of the Western Blot was performed by densitometry, using ImageJ software. 
doi:1 0.1 371/journal.pgen.1 0041 55.g001 
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Figure 2. Features of SREBP1 binding sites. (A) Distribution of peak lengths in the four clusters of SREBP1 binding sites shown in figure 1 A. (B) 
Distribution of the distance from the closest annotated TSS of SREBP1 binding peaks belonging to the four clusters. Cluster A is enriched in 
sequences closer to a TSS, whereas the profile of cluster B, C and D is overlapping with that of 1000 randomly selected sequences (black line). (C) 
Distribution of the amplitudes of SREBP1 binding oscillation. (D) Overrepresented motifs within SREBP1 binding sites belonging to cluster A were 
found using IVIEIVIE [21]. DNA sequences in the window under each SREBP1 peak were used for the sequence analysis. (E) The 236 SREBP1 binding 
sites (distributed in 223 genes) belonging to cluster A were associated to the closest genes. Annotations are available for 219 of these 224 genes and 
the Venn diagram shows the overlap between the presence of SREBP1 and the presence of a nearby a binding site for SP1, NFY and/or HNF4. The 
most relevant functional pathways that were enriched in the different sets of target genes are indicated. 
doi:1 0.1 371/journal.pgen.1 0041 55.g002 



transcript. We used DAVID [22,23] to identify clusters of genes 
enriched with functional annotations. As expected, we identified 
lipid biosynthetic processes and fatty acid metabolism as the most 
prominent pathways controlled by SREBPl (Table S3). In 
addition, we found a significant enrichment of genes involved in 
carbohydrate metabolism, in the response to nutrient levels, in 
mitochondrial and endoplasmic reticulum functions and in 
coenzyme metabolism. In a previous genome-wide study per- 
formed in human HepG2 cells, it was shown that unique 
combinations of SREBPl, SPl and NFY target distinct functional 
pathways [17]. Since we found a good enrichment within the 
SREBP 1 binding sites of the consensus motifs for NFY and SP 1 , 



but also HNF4, we explored whether a network among these three 
transcription factors could be highlighted in mouse liver. As shown 
in Table 1, genes involved in lipid biosynthesis and in the 
regulation of fatty acid and steroid metabolism were highly 
represented in all categories. In some cases, however, one 
biological function was targeted by a unique combination of 
regulators. For example, the biosynthesis of coenzymes was 
selectively represented within the genes bearing only the site for 
SPl, whereas both SPl and HNF4 motifs were present in genes 
involved in apoptosis. Likewise, a combination of HNF4 and NFY 
motifs marked most of the genes involved in immunological 
processes. Finally, pathways related to carbohydrate metabolism 
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Table 1. Functional annotation clustering of putative SREBPl 
targets using DAVID tools. 





Binding Motif 


n 


Functional cluster 


GO Terms 


P-Value 


SPl 


65(63) 


sterol metabolism 


GO:0016125 


6.50E-07 






coenzyme biosynthesis 


GO:0009108 


1 .40E-03 






metal ion binding 


GO:0046872 


6.30E-04 






lipid biosynthesis 


GO:0008610 


4.20E-04 






magnesium ion binding 


GO:0000287 


4.20E-05 






cell fractions 


GO:0000267 


0.029 


NFY 


21(20) 


lipid biosynthesis 


GO:0008610 


2.70E-06 


HNF4 


22(22) 


No clusters 






SPl+NFY 


32(32) 


nucleotide binding 


GO:0005524 


2.70E-03 






sterol metabolism 


GO:0016126 


9.60E-04 






inflammatory response 


GO:0006954 


0.04 


SP1+HNF4 


18(18) 


cholesterol metabolism 


GO:0008203 


2.80E-05 






regulation of apoptosis 


GO:0042981 


0.033 


NFY+HNF4 


15(14) 


antigen processing.... 


GO:0002474 


6.60E-05 


SP1+NFY+HNF4 


4(3) 


No clusters 






Only SREBPl 


65(63) 


glycerolipid metabolism 


GO:0046486 


9.40E-03 






fatty acid metabolism 


GO:0006631 


0.024 






glucose metabolism 


GO:0006006 


0.041 






mitochondria 


GO:0005739 


5.00E-03 



Enriched GO categories were identified in four distinct sets of SREBPl target 
genes exhibiting a different combination of binding sites for SPl, NFY and/or 
HNF4. In total, 219 out of 223 SREBPl putative target genes have a functional 
annotation (the number of annotated genes for each set is indicated in 
parentheses}. The analysis using DAVID groups the GO categories in functional 
related clusters. For each enriched cellular process, only the most significant 
associated GO term is shown in the table, with the corresponding Modified 
Fisher Exact P-value. The complete list of all GO terms enriched in each 
functional cluster for all the groups is available in Table S4. 
doi:l 0.1 371 /journal.pgen.l 0041 55.t001 

and mitochondria were particularly enriched in genes without 
NFY or SPl motifs, suggesting that SREBPl may cooperate with 
other regulators at the promoter of these genes (the complete 
functional annotation clustering is reported in Table S4). Our 
results suggest that in mouse liver, in physiological conditions, the 
network SREBP1-SP1-NFY-HNF4 may be important in order to 
determine the functional effect of SREBPl binding (Figure 2E). 

RNA polymerase II recruitment to SREBPl target genes is 
not always synchronized with SREBPl binding 

In Figure 1, we showed that SREBPl binds to target sites 
belonging to cluster A with a sharp phase between ZT15 and 
ZT 1 7 . To investigate the functional effects of SREBP 1 binding on 
gene transcription, we checked in our previously reported data set 
[24] the 24 hours profile of RNA polymerase II (Pol II) 
recruitment in the proximity of SREBPl target genes. Important- 
ly, most of the SREBP 1 binding sites belonging to cluster A were 
co-occupied by Pol II (Figure lA), further supporting the 
functional relevance of these regions. We next evaluated Pol II 
binding to the promoter and in the gene body of all putative 
SREBPl target genes. In parallel, we measured mRNA levels of 
the same genes by microarray analysis (Table S5). More than 85% 
of SREBPl target genes show an expression level above the 
median expression level of aU the transcripts, suggesting that they 
are transcribed. Our analyses revealed three clusters of target 



genes, that we called Al, A2 and A3, with distinct temporal profile 
of transcription and expression (Figure 3A and Figure S2). In 
cluster Al, the peak of Pol II binding was concomitant, or even 
slightly earlier than SREBPl binding. In contrast, for genes 
belonging to cluster A2, Pol II association to both promoter and 
gene body strictly followed SREBPl binding. Lastly, Pol II 
recruitment to the genes of the A3 group was shifted by about + 
8 h with respect to SREBPl. For all clusters, the temporal profile 
of gene expression was consistent with the dynamics of Pol II 
association. The distribution of aU the expression phases obtained 
for the genes belonging to the three groups confirmed that 
SREBPl target genes are expressed in different moments of the 
day, in spite of the concomitant binding of the transcription factor 
(Figure 3B). This observation suggests that other factors participate 
in the regulation of the various SREBP 1 target genes in order to 
assure their appropriate expression timing. Interestingly, genes 
that were mainly expressed during the fed state (clusters Al and 
A2) were functionally enriched in the regulation of lipid and 
coenzyme biosynthetic processes, as well as in the response to 
hormones, such as insulin. In contrast, SREBPl target genes 
involved in mitochondrial oxidation and apoptosis were enriched 
during the fasting period (Table 2). Thus, the promoter specific 
events that determine the different temporal expression profile of 
SREBPl target genes contribute to define the set of cellular 
functions that are active at a given time. 

A functional core molecular clock is necessary for the 
timing of SREBPl target gene expression 

To understand the molecular mechanism underlying the 
different temporal expression of SREBPl target genes, we first 
explored the possible involvement of the network SREBP 1 -SP 1 - 
NFY-HNF4 in determining the functional effect of the binding of 
SREBPl to its targets. To check for the presence of a pattern 
characterizing the three groups of SREBP 1 target genes identified 
earlier (see Figure 3), we evaluated the presence of different 
combinations of SPl and NFY motifs and their orientation with 
respect to the SREBPl binding sites (data not shown). However, 
we could not establish any significant correlation. In contrast, we 
found that HNF4 motifs were significantly overrepresented (P- 
value<0.02) in the regions under SREBPl peaks of the genes 
expressed during the fasting period (cluster A3), compared to the 
other clusters (Table S6). The actual recruitment of HNF4 to these 
putative binding sites was assessed by ChIP on randomly selected 
SREBPl Responsive Elements (SREs) (Figure 3C). Besides HNF4, 
other transcription factors, such as the cAMP response element- 
binding protein (CREB) or Forkhead box proteins O (Foxo), are 
important players in the hepatic metabolic regulation upon fasting 
[25,26]. However, their known consensus motifs were not found in 
the proximity of cluster A3 SREBPl peaks. These observations 
strongly suggest a specific cross-talk between HNF4 and SREBPl 
in the regulation of these genes. To further investigate which 
control processes dictate the distribution of SREBPl target gene 
expression along the day, we then considered the possible role of 
the circadian rhythm in this regulation. To test this hypothesis, it 
was necessary to uncouple the circadian rhythm from the response 
to nutrients. Thus, we fed mice lacking BMALl [Bmall~ ^ ~) only 
during the darkness period for one week before collecting liver 
samples every four hours. Upon this experimental conditions, the 
circadian clock was completely disrupted in Bmall ' mice, as 
demonstrated by the flattened expression of key core and output 
components of the clock, such as Clockl, Cryl, Cr)i2, D site albumin 
promoter binding protein [Dbp), Rev-Erba and Kruppel-like factor 
10 [KlflO) (Figure 4B). Body weight, daily food intake and glycemia 
were unchanged in Bmall~'~ mice (Figure S3). Importantiy, the 
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Figure 3. Pol II recruitment on SREBP1 target genes is not always synchronized with SREBP1 binding. (A) The heat-maps represent the 
recruitment of Pol II to the promoter (left) and to the gene body (middle) of SREBP1 putative target genes along the day, as assessed in our previous 
ChlP-seq data set [24]. In parallel, we evaluated hepatic gene expression by microarray analysis in the same samples (right). Hierarchical clustering 
was done applying a Pearson correlation scores to the data describing Pol II recruitment on the promoters of SREBP1 target genes (left). Three major 
clusters of genes displaying a different temporal binding profile of Pol II were identified (A1, A2, and A3). The genes are ordered in the three heat- 
maps according to this clustering. (B) Gene expression data from microarray analysis were fit to a cosine function to estimate the phase of expression 
(peak time of the fit) of SREBP1 target genes. The graph shows the smoothing of phase distributions of the genes belonging to the three clusters 
(green line for A1, blue line for A2, magenta line for A3). Only genes with a P-value<0.05 are plotted. Dotted lines define three time intervals 
containing the most recurrent phases associated to the genes belonging to the clusters A1, A2 and A3. (C) HNF4 binding was tested on randomly 
selected SRE identified in our SREBP1 ChlP-seq. Gnatl and AnksAb SREs belong to cluster A3 and contain a HNF4 putative binding sites. In contrast, 
Ldir and Insigl SREs do not contain in their sequence a HNF4 motif and belong to cluster A2. NEG and POS were used as negative and positive control 
loci and correspond to two regions of Cyp7al promoter, localized at -1500 and -150 from the TSS, respectively [69]. The graph shows the mean ± 
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SEM of three independent experiments. * indicates P-value<0.01 vs. NEG. Statistical analysis was performed by one-way ANOVA followed by 
Bonferroni post-test. Primer sequences are listed in Table S8. 
doi:1 0.1 371 /journal.pgen.1 0041 55.g003 



imposed rhythmic food intake restored an oscillatory nutrient 
response, as shown by the levels of circulating insulin that 
were comparable to the wild type (Figure 4A). Accordingly, in 
Bmall~'~ mice SREBPl translocated to the nucleus from ZT18 
onwards (Figure 4D and 4E) and its expression was stiU cycling, 
although the phase was delayed by about 6 h compared to the wild 
type (Figure 4C). We next checked the dynamics of SREBPl 
binding to a panel of the targets previously identified in wild type 
mice and we found that SREBPl was recruited to all the tested 
sites in an oscillatory way, but with an average phase shift of about 
4 hours (figure 4F). Finally, to assess the impact of the circadian 
oscillator impairment on SREBP 1 -driven transcription, we glob- 
ally evaluated the expression of SREBP I target genes in Bmall ' 
mice (Table S5). Most SREBPl target genes were still scored as 
oscillating in Bmall ' upon temporal restricted feeding (57% 
have a P<0.05). However, the heatmap rendering of their 
expression patterns (Figure 5A) revealed a temporal profile that 
was perturbed in Bmall~^~ compared to WT mice, as most of the 
genes now had a maximum expression at ZT18, coinciding with 
the binding of the transcription factor (Figure 5C). Accordingly, 
the expression phases of the genes belonging to the clusters Al and 
A3 identified earlier (Figure 3) were now largely concomitant with 
those of genes belonging to cluster A2, therefore mosdy grouped 
between ZT14 and ZT24 (Figure 5B and Figure S4). This phase 
shift was not due to a selective decrease of the number of cycling 
genes in clusters Al and A3, as the percentage of the significandy 
oscillating genes was comparable in the three clusters in WT and 
Bmair'~ mice (cluster Al: 84% in WT vs 94% in Bmair'~; 

Table 2. Functional annotation clustering of putative SREBPl 
targets with a different temporal expression profile. 



PHASE 


n 


Functional cluster 


GO term 


P-value 


ZT13-ZT19 


42 


cholesterol metabolism 


GO:0008203 


0.01 






amino acid binding 


GO:0016597 


2.00E-03 






response to starvation 


GO:0042594 


1 .90E-03 






response to hormones 


GO:0009725 


6.10E-03 






glycerolipid metabolism 


GO:0046486 


3.10E-03 






cell migration 


GO:0016477 


0.017 


ZT19-ZT28.4 


68 


sterol metabolism 


GO:0016126 


1 .60E-09 






magnesium ion binding 


GO:0000287 


1 .OOE-03 






coenzyme biosynthesis 


GO:0006732 


0.019 






nucleotide binding 


GO:0032555 


0.041 






endoplasmic reticulum 


GO:0005783 


0.041 


ZT4.4-ZT13 


39 


mitochondria 


GO:0005739 


3.20E-03 






regulation of apoptosis 


GO:0042981 


6.70E-03 






regulation of myeloid cell 
differentiation 


GO:0045637 


2.30E-03 



Enriched GO categories were identified in three distinct sets of rhythmic SREBPl 
target genes (P<0.05), based of their phase of expression. To define the three 
intervals of time we calculated the shortest time range containing the phases of 
at least 50% of the genes belonging to clusters Al, A2 or A3. For each set, the 
total number of genes is indicated and the number of genes with annotation is 
indicated in parentheses. For each functional cluster, only the most significant 
associated GO term is shown in the table, with the corresponding Modified 
Fisher Exact P-value. 
doi:l 0.1 371/journal.pgen.1 0041 55.t002 



clmter A2: 71% in WT vs 66% in Bmall ' ; cluster A3: 62% in 
WT vs 64% in Bmair'~ with P<0.05). 

To explore how the core clock components participate to this 
regulation we checked in published data sets whether key 
transcription factors such as BMALl, CLOCKl, CRYl, CRY2, 
PERI, PER2, NPAS2 and REV-ERBs can diflferentiaUy bind to 
the promoters of the three clusters of SREBPl target genes [27- 
29]. Interestingly, we found that most SREBPl peaks (148 out of 
236, «63%) have an overlapping REV-ERBa and/or REV- 
ERBP peak. In addition, a REV-ERBa/p binding site was 
detected also in another 1 7 % of the promoters of SREBP 1 target 
genes, but in a non-overlapping position. This strong occupancy of 
SREBPl targets by REV-ERBs is consistent with the previously 
reported involvement of REV-ERBa in the regulation of lipid 
metabolic genes [30], and suggests the existence, in mouse liver, of 
a SREBPl -REV-ERBs network in physiological conditions. The 
frequency of REV-ERBs recruitment was comparable in clusters 
Al, A2 and A3 (data not shown), thus arguing against the possible 
role of these nuclear receptors in determining the distinct phase of 
expression of these genes. However, due to the presence of REV- 
ERB binding sites in many SREBPl target genes, the flattened 
REV-ERB expression observed in Bmall ' may perturb, at least 
in part, die phase of several SREBPl target genes. Indeed, in WT 
mice, the temporal expression profile of SREBP 1 and REV-ERBs 
is very dififerent (the phases of expression are ZT15 and ZT8, 
respectively [16]), and these factors are not expected to compete 
for binding at the same time to the same genes. The other 
transcription factors tested were recruited to a lesser extent on 
SREBP 1 target gene promoters and for none of them we observed 
a significant enrichment in clusters Al, A2 or A3 (data not shown). 

Taken together, our results confirm that SREBPl activity is 
strongly dictated by the rhythmicity of nutrient intake. In addition, 
our observations indicate that a functional circadian core clock is 
necessary to assure the correct temporal expression profile of 
SREBP 1 target genes and suggest a role for HNF4 in dictating the 
phase of expression of genes whose mRNA levels peak when 
SREBP 1 binding is low. Further studies will aim at understanding 
whether and how the lack of circadian rhythm perturbs HNF4 
activity. 

Discussion 

SREBP 1 is a highly circadian transcription factor whose activity 
is strongly regulated by nutrient availability through the insulin 
signaling pathway. In mouse liver SREBPl expression displays a 
daily rhythm with a peak in the nocturnal feeding period under 
standard housing condition of mice [16,18-20]. In this study we 
evaluate the dynamics of SREBPl recruitment to DNA by 
determining its genome wide cis-acting targets (cistrome) in the 
liver along an entire day. SREBPl binds to 448 sites with an 
oscillatory profile that is temporally coherent with the phase of its 
maximal expression. Within SREBPl binding sites, four distinct 
groups are clearly distinguishable. The first set (cluster A) contains 
peaks that are likely the more relevant in the transcriptional 
regulation mediated by SREBPl as they are the closest to TSS and 
they are bound more rhythmically by SREBPl. Importantiy, in 
more than 60% of these sites we identified the direct repeat 5'- 
ATCACCCCAC-3' that was described as the Sterol regulatory 
proteins Responsive Elements (SRE) in several promoters, such as 
the human LDL receptor promoter [13,31-33]. This direct repeat 
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Figure 4. SREBP1 binding is rhythmic in Bmah '^ upon time-restricted feeding. (A) Bmall^^^ (red line) and control mice (black line) were 
fed only during the night for one week before the sacrifice. Plasma insulin levels were measured at the indicated time points (n = 3-6). (B) mRNA 
levels of key genes of the cellular molecular clock were measured by qPCR in Bmall^^^ and control mice (n = 5). (C) Hepatic expression of Srebpl c 
was evaluated by qPCR in Bmall^^^ and control mice (n = 5). Data are normalized using 36b4 and /?ps9 as housekeeping genes. (D) Representative 
western blot analysis of the nuclear SREBP1 in hepatic nuclear extracts from Bmall^^^ mice. Lamin A was used as loading control. Each sample is a 
pool of 5 livers. (E) Western Blot quantification was performed by densitometry, using ImageJ software. Each point represents the mean ± SEIVI of the 
quantification of three analyses performed in three independent sets of Bmall^^^ mice (n>3 per each time point). (F) ChIP of SREBP1 was performed 
in livers of Bmall^^^ and control mice at the indicated time points. SREBP1 binding was tested on 6 loci (SRE) identified by ChlP-seq in the proximity 
of the indicated genes. Aacs and Srebpl c belong to cluster A1, Insigl and Clcn6 belong to cluster A2, whereas Obfc2a and Slc25a25 belong to cluster 
A3. Primer sequences used for qPCR analyses are available in Tables S7 and S8. 
doi:1 0.1 371 /journal.pgen.1 0041 55.g004 
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Figure 5. Phase of expression of rhythmic SREBP1 target genes is perturbed in Bmall upon time-restricted feeding. (A) The heat- 
maps represent the hepatic gene expression of putative SREBP1 target genes In wild type (left) and Bmall^^^ mice (right), after one week of time- 
restricted feeding, as assessed by microarray analysis. The order of the genes Is based on the clustering shown In figure 3 and Is maintained In the 
two heat-maps. (B) Validation of hepatic mRNA level variation of a panel of the Indicated SREBP1 putative target genes. (C) Gene expression data 
from microarray analysis In Bmall^^^ mice were fit to a cosine function to estimate the phase of expression (peak time of the fit) of SREBPl target 
genes. The graph shows the smoothing of phase distributions of the genes belonging to the three clusters (green line for A1, blue line for A2, 
magenta line for A3). Only genes showing a P value<0.05 in Bmall^^^ mice are plotted. 
dol:1 0.1 371 /journal.pgen.1 0041 55.g005 

variant of the canonical E-box inverted repeats 5'-CAnnTG-3' an arginine in all the other bHLH-LZ proteins and that is critical 
was shown to be specifically recognized by SREBP proteins due to for high affinity contacts with the SRE [34,35]. Furthermore, our 
the presence of a tyrosine residue in a position that corresponds to ChlP-seq results highlighted the presence of predicted binding 
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sites for SPl, NFY and HNF4 in 60%, 30% and 25% of cluster A 
sites, respectively. Since SREBPlc, the major SREBPl isoform in 
the liver, is a weak transcriptional activator, this observation is 
consistent with earlier studies demonstrating that the transcription 
factors SPl, NFY and CREB cooperate to regulate different 
SREBPl -responsive promoters [36-39]. In 2009, Seo et al have 
published a list of liver SREBPl target genes obtained from a 
genome-wide study of mice subjected to 24 hours fasting followed 
by 12 hours refeeding with high carbohydrate diet [40], thus 
creating a condition where a very high SREBPlc activity is 
expected. In this study, a funrtional variant of the direct repeat 
SRE (5'ACTACANNTCCC-3') was identified as a preferred site 
for SREBPl binding, and no enrichment of the predicted NFY 
binding site was identified. The difference between this data set 
and ours can most likely be attributed to the difference in the 
specific experimental conditions, acute challenge on the one hand 
and physiological condition on the other hand (present study). 
Consistent with our study, both SPl and NFY proteins were 
recruited on more than 30% of SREBPl target genes in a genome- 
wide analysis of SREBPl binding in HepG2 cell line [17]. In 
addition to SP 1 and NFY, here we identify HNF4 as an important 
player of the interconnected regulatory circuit that may assure the 
specific regulation of SREBPl target genes with distinct functions. 
Interestingly, consensus motifs for SPl, NFY and HNF4 were 
found to be overrepresented in the promoters of cycling genes in 
the liver [41,42]. Since =70% of SREBPl targets show a circadian 
gene expression, our results are in line with these bioinformatics 
predictions, supporting the involvement of these transcription 
factors in the complex transcriptional regulation of circadian 
rhythm in liver. The second set of SREBP 1 target sites falls in the 
three clusters B, C and D. These peaks are not enriched in regions 
proximal to TSSs for mapped genes, nor in predicted motifs for 
known transcription factors. Furthermore the temporal profile of 
SREBP 1 binding to these sites is flattened compared to the first set 
of SREBPl target sites (cluster A). Additional studies are required 
to understand whether these peaks have a functional role, or 
whether they are bound in a secondary manner by SREBPl due to 
the formation of DNA loops. 

In the liver the expression of several known SREBPlc target 
genes is decreased in fasted mice, when the levels of SREBPl arc- 
very low and increased upon refeeding, when both SREBPl 
expression and nuclear translocation are induced [40,43]. 
Accordingly, our analysis of Pol II recruitment on SREBPl 
putative target genes, coupled with the measurement of their 
mRNA levels, revealed a maximum of transcription and expres- 
sion during the fed state, namely between ZT 1 2 and ZT24, for the 
majority of these genes. This is consistent with the binding of 
SREBPl to DNA, which is higher at this time of the day. 
Nevertheless, a large set of SREBPl target genes (cluster A3) 
displays a temporal expression profile strongly shifted with respect 
to SREBPl recruitment. However, the different phase of 
expression observed in these genes is coherent with the dynamics 
of Pol II association to their promoter and gene body, thus arguing 
against a major involvement of post-transcriptional mechanisms in 
this delay and suggesting the existence of promoter specific events 
that determine the different temporal expression profile of 
SREBPl target genes. In particular, the low level of expression 
of genes belonging to cluster A3 when SREBP 1 is bound raises the 
question whether SREBPl itself, or through interaction with 
coregulatory proteins, can act as transcriptional repressor for these 
genes. Interestingly, it was proposed that SREBPl may act as 
negative regulator of the cytosolic phosphoenolpyruvate carbox- 
ykinase (Pck-1) gene by impairing the recruitment of the 
transcriptional coactivator Peroxisome Proliferator-Activated 



Receptors y Coactivator -1 (PGC-1) on HNF4a [44]. To (-xplain 
the negative effect of SREBPl on this gene, a second mechanism 
was put forward by Chakravarty and colleagues, suggesting an 
interference between the binding of SREBPl and SPl, due to the 
orientation on the opposite DNA strands of the two binding sites 
[45]. Although we find an enrichment of motifs for NFY and SPl 
in the close proximity of SREBPl peaks from the first set of target 
sites (clusters Al, A2 and A3), we do not observe a difierent 
presence and/or orientation of these sites with respect to SREBPl 
peaks, among these three groups of target sites. Conversely, the 
high frequency of the HNF4 motif in the cluster A3 suggests that 
the cross-talk between HNF4 and SREBPl may be a general 
mechanism through which SREBPl negatively affects the 
transcription of a sub-set of its target genes. In agreement with 
this hypothesis, among the targets of SREBPl expressed upon 
fasting we detect the Peroxisome Proliferator-Activated Receptor 
a iPpara) gene, that was shown to be crucially regulated by HNF4 
[46]. 

In recent years, growing evidences have highlighted the impact 
of circadian gene networks on nutrient balance and, on the other 
hand, the regulation of the circadian clock by metabolism and 
food consumption [47,48]. Thus, circadian clock and metabolism 
converge in numerous ways to control the activity of a number of 
transcription factors that are essential for maintaining metabohc 
homeostasis, although the exact contribution of each input 
remains to be deciphered. Several studies demonstrated that upon 
restricted feeding (RF), namely when time and duration of food 
availability is limited in time, mice adjust to the feeding period 
within a few days, they display food anticipatory behavior and 
consume their daily food intake during that limited time [49-52]. 
This feeding regimen drives rhythms in arrhythmic and clock 
mutant mice or in animals with SNC ablations, thus uncouphng 
the circadian clock, .synchronized by SCN, from the periphery 
[3,4] . Many physiological activities that are normally dictated by 
the SCN master clock, such as hepatic P450 activity, body 
temperature, locomotor activity, and heart rate, are restored by 
RF. In the liver of cryl ' ;ciy2 ' mice, RF restores the 
oscillatory circadian expression profile of a number of "feeding 
driven" transcripts, although with a small delay in their phase of 
expression, showing that the circadian (:lot:k anticipates changes in 
the feeding state and accelerates the transcriptional response to an 
acute activation or repression by feeding [4]. This is consistent 
with our observation that in Bmall ' mice, rhythmic SREBPl 
expression and activity, that are drastically flattened when mice 
are fed ad libitum (data not shown), are reinstated upon RF, 
although with a deferred phase. Notably, growth and metabolic 
defects that were reported in Bmall~'~ mice either at older age or 
under different feeding regimens [53-58], are negligible in the 
experimental conditions adopted in our study, suggesting that the 
role of the circadian clock in the regulation of SREBPl could be 
evaluated in the absence of major confounding pathologies. As an 
example, Bmall ' mice at 8-10 month of age have an impaired 
insuhn release due to the absence of a functional clock in 
pancreatic beta-cells [59]. However, in our case the delayed 
SREBPl activation cannot be attributed to a reduced insulin 
release, since we detect normal glucose and insulin levels in 
Bmall~'~ mice at 3 month of age upon RF. Conversely, the 
expression of REV-ERBot, that is directly regulated by BMALl, is 
constantiy downregulated. This event leads to the derepression of 
the Insulin Induced 2 (InsigZ) gene, encoding a trans-membrane 
protein that sequesters SREBP proteins to the endoplasmic 
reticulum membranes, thus interfering with the proteolytic 
activation of SREBPs, in agreement to what was shown in 
REV-ERBa-/- mice [16]. SREBPl activity in the nucleus 
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reflects also the rate of its proteosomal degradation after DNA 
binding [60,61], a process that is strongly sensitive to the insulin- 
mediated inactivation of the glycogen synthase kinase 3 (Gsk3) 
[62]. In Bmall ' mice the phase of SREBPl recruitment to 
DNA is shifted, but we do not observe a longer SREBPl 
accumulation on its targets, suggesting that the absence of a 
functional clock is not significantly altering its degradation process. 
In conclusion, our results show that besides the nutrient-driven 
regulation of SREBPl nuclear accumulation, a second layer of 
modulation of SREBPl transcriptional a( ti\'ity exists and is 
strongly dependent from the circadian core clock. This system 
allows to fine tune the expression timing of SREBPl target genes, 
thus helping to temporally separate the different physiological 
processes in which these genes are involved. Thus, SREBPl is 
situated at the interface of the circadian and the metabolic 
regulation and its study promises to shed light on the emerging 
association between diabetes, obesity, sleep, and circadian timing. 

Materials and Methods 

Animals 

AU animal experiments and procedures were approved by the 
Swiss Veterinarv- Office (authorisation VD-14.53.4). C57BL/6 
male were purchased from Charles River. Bmall ' mice were a 
kind gift from Dr. Frederic Gachon and were generated as 
previously described [63,64]. 12-14 week old (at time of sacrifice), 
mice were housed in a 12 h light/ 12 h dark (LD) regimen for 2 
weeks with food and water freely available during night and day. 
They were then phase-entrained to a 1 2 hr/ 1 2 hr LD regimen 
with food access between ZT12 and ZT24 for 7 days 
(ZT = Zeitgeber time; ZTO is defined as the time when the lights 
are turned on and ZT12 as the time when lights are turned oil). At 
each ZT2, ZT06, ZTIO, ZT14, ZT18 and ZT22 three to five mice 
were anesthetized with isoflurane and decapitated. Mice were 
killed under dim red light at ZTs during the dark phase. The livers 
were perfused with 2 ml of PBS through the spleen and 
immediately collected. A small piece of liver tissue was snap- 
frozen in liquid nitrogen. The remaining liver tissue was 
immediately homogenized in PBS containing 1% formaldehyde 
for chromatin preparation. 

Chromatin immunoprecipitation (ChIP) 

Perfused livers were processed for chromatin preparation as 
previously described [65]. The chromatin samples from the mice 
of the same ZT were then pooled, frozen in liquid nitrogen and 
stored at — 80°C. The following antibodies were used: anti-RPB2 
(Santa Cruz Biotechnology, H-201), anti-SREBPl (Santa Cruz 
Biotechnology, H-160), anti-HNF4 (Santa Cruz Biotechnology, C- 
19). Chromatin was subjected to immunoprecipitation of Pol II as 
described [24] . For SREBP 1 , the samples were diluted ten times in 
"sonication buffer" containing 50 mM HEPES (pH 7.9), 140 mM 
NaCl, 1 mM EDTA, 1% Triton X-100, 0.1% Na-deoxycholate, 
0.1% SDS and proteinase inhibitors (Roche). 1 ml of diluted 
chromatin was immunoprecipitated with 10 ng of antibody as 
described [66]. Briefly, the immune complexes were collected by 
adsorption to ten |al of protein-A-Sepharose (25% slurry in 
sonication buffer), pre-blocked with 10 |J.g/ml of salmon sperm 
DNA and BSA at 4°C overnight. The beads were washed twice 
with "sonication buffer", twice with sonication buffer containing 
500 mM NaCl, twice witii 20 mM Tris, pH 8.0, 1 mM EDTA, 
250 mM LiCl, 0.5% NP-40, 0.5% Na-deoxycholate and twice 
with TE buffer. The immunocomplexes were eluted with 50 mM 
Tris, pH 8.0, 1 mM EDTA and 1% SDS at 65°C for 10 min., 
adjusted to 200 mM NaCl and incubated at 65 °C overnight to 



reverse the cross-links. After successive treatments with 10 ng/ml 
Rnase A and 20 |Xg/ml proteinase-K, the samples were extracted 
with NucleoSpin Kit (Macherey- Nagel). The DNA concentration 
was determined by fluorometry on the Qubit .system (Invitrogen). 
10-12 ng DNA were used for the preparation of the library. 
Libraries for ultra-high throughput sequencing were prepared with 
the ChlP-Seq DNA sample kit (lUumina) as recommended by the 
manufacturer. 

RNA isolation and analysis 

About 100 mg of snap-frozen liver tissue were used for RNA 
preparation with the TRIzol reagent (Invitrogen) followed by 
purification with miRNeasy Mini Kit (Qiagen), according to 
manufacturer's instructions. For microarray analysis 500 ng of 
total RNA from each liver sample at the same time point were 
pooled and analyzed on Mouse Gene LOST arrays according to 
the manufacturer's instructions (Aflymetrix). All statistical analyses 
were performed with the statistical language R and various 
Bioconductor packages (http://www.Bioconductor.org). Normal- 
ized expression signals were calculated from Affymetrix CEL files 
using RMA normalization method. For quantitative RT-PCR 
analysis, the retrotranscription has been done using iScript cDNA 
synthesis kit (Bio-Rad, Laboratories, Hercules, CA) and following 
the manufacturer's instructions. The primers sequences are shown 
in Tables S7 and SB. Real-time monitoring of PCR amplification 
of cDNA was performed using the FastStart Universal SYBR 
Green Master (Roche Applied Science, Indianapolis, IN) in an 
ABI Prism 7900 Sequence Detection System (Life Technologies, 
Carlsbad, CA). The PCR arbitrary units of each gene were defined 
as the mRNA levels normalized to the 36h4 and the Rps9 
expression level in each sample using the qBase Software. 

Western blotting 

Nuclear extracts were prepared by the NUN procedure as 
described previously [67], and Western blotting was performed 
according to standard protocols using the antibody for SREBPl 
indicated above. U2AF and Lamin A were used as loading control 
(anti-U2AF and anti-Lamin A were from Sigma-Aldrich). 

Blood biochemistry 

At sacrifice, blood was taken for determination of biochemical 
parameters and circulating hormones. Insulin levels were deter- 
mined with ELISA kit from Mercodia, Uppsala, Sweden, 
foUowing manufacturer's instructions. 

ChlP-seq data analysis and read mapping 

At eat:h time point, DNA sequenced reads were mapped to the 
mouse genome (Mus musculus NCBI m37 genome assembly 
(mm9; July 2007)) using Bowtie [68] with three mismatches and at 
most five hits aUowed on the genome. When computing genomic 
read densities, each alignment contributed 1 /(total number of hits) 
to the local density. If several reads in the same library- mapped at 
the same genomic position and on the same strand (redundant 
tags), we kept only one read for the rest of the analysis. The total 
numbers of reads per time point are given in Table SI. 

Strand shifting 

The mapped reads were shifted to account for the length of the 
inserts based on the average fragment size, namely 190 184 204 
199 202 208 and 176, for each of the seven libraries. The fragment 
size was divided by two and half of the read length was subtracted 
resulting in a shift of 55, 52, 62, 60, 61, 64, and 48 nucleotides 
from ZT02 to ZT26. 
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Peak identification and refinement 

We fragmented the genome into 500 nucleotide blocks and 
collected the counts within each block for SREBPl as well as for 

input (experiments. We kept all blocks for which we had a signal 
equivalent to 40 tags in at least one time point of tlie SREBPl 
experiment. We log transformed the data after adding 1 pseudo 
count and quantUe normalized both SREBPl and input experi- 
ments. We selected blocks with a log2 signal SREBPl /input 
greater than 2, i.e. at least a four-fold enrichment of the SREBPl 
signal in comparison to the input signal in at least one time point. 
We repeated this procedure by shifting the block definition of half 
the length of the blocks and merged the overlapping blocks that 
passed these criteria. 

To define proper "peaks" in these wide regions, we look for 
shorter regions accounting for most of the counts. For this, we 
repeatedly consider the two borders along 50 nucleotides, and 
discard the one with less read counts if we keep 75% of the total 
reads in the remaining region. This operation is repeated with 
shorter borders until no further refinement is possible. 

Motif search 

We used the MEME suite [2 1] to identify enriched motifs in the 
sequences corresponding to the refined peaks. Wt; first performed 
several motif discovery, on all peaks, and on the subset associated 
or not with Pol II. We searched for 15 motifs between 6 and 10 
nucleotides long. The discovered motifs have been associated to 
known transcription factors (from the TRANSFAC database) with 
STAMP. To retrieve these motifs in the subsets where they have 
not been discovered, we searched for them in all sequences (using 
FIMO). To assess the relevance of the number of observed motifs 
in our dataset, we counted the occurrence of the same motifs in 
random sequences. These sequences are selected in the proximity 
of the TSS of genes expressed in our samples (among the top 10% 
in the microarray data). Their size and distance to TSS are in the 
same range as that of our SREBPl peaks. 

Cosine fits 

Before fitting a cosine function to estimate the amplitude and the 
phase of the oscillation in a 24 hour period, counts in the refined 
peaks were quantified and normalized according to the total 
number of non redundant mapped reads for each given library. 

We used the function x(t) = b0-l-bl*cos(b3-l-27l*t/24) to perform 
a least squared fitting of temporal profiles. The parameter bO 
represents the mean signal, b 1 the amplitude of the oscillation, and 
b3*24/2Ti the phase. These parameters were estimated by 
nonlinear least-squares using the Gauss-Newton algorithm. 

For the microarray temporal profile analysis we used the 
function 

x(t) = bO-l-blcos(b3-l-27t * t/24) -|-b4BE, where b4 represent 
the batch effect 

while when we compared WT and Bmall ' mice we used the 
function 

x(t) =bO-l-bOgtGT-|- (bl -hblg,GT)* 

cos[(b3-hb3gtGT)27c * t/24] -hb4BEl -hb5BE2, 

where b4 and b5 represent the different batch elfects, GT is a 
dummy variable that indicates the mouse genotype (WT or KO) 
and bOgt, blp, b3g, the associated coefiicients. 

Data availability 

lUumina sequencing data for the ChlP-seq are available at 
GEO as the GSE48375. Additional processed data and visualiza- 
tion tools are provided at http:/ /cyclix.vital-it.ch. 



Supporting Information 

Figure SI VaUdation of SREBPl antibody for ChIP experi- 
ments. (A) C57BL/6 mice were fasted 24 hours and re-fed for 
12 hours before liver collection for chromatin preparation. 
SREBPl binding was tested on two positi\'c control loci, the Fatty 
acid synthase {Fasrij and Low density lipoprotein receptor [Ldlr) promoters. 
Negl and Neg2 were used as negative control loci, and correspond 
to a site specifically recognized by Pol III and to a region between 
a exon 6 and intron 6-7 of the glyceraldehyde 3-phosphate dehydrogenase 
(Gapdh) gene, respectively. Fold enrichments relative to the 
negative controls are greater than 40-fold for Fasn and about 10- 
fold for Ldlr. (B) In samples processed in absence of antibody none 
of the tested sequences was significantly enriched. Primer 
sequences are listed in Table S8. 
(PDF) 

Figure S2 Temporal relationship of SREBP 1 and Pol II profiles 
and mRNA accumulation of SREBPl putative target genes. Phase 
histograms for SREBPl target genes belonging to clusters Al, A2 
and A3. While SREBPl binding phases are sharply concentrated 
between ZT14 and ZT18 for all genes, the temporal distribution 
of Pol II recruitment and mRNA expression is different in the 
three clusters. The phases are plotted only for genes with an 
amplitude P-value<0.05 for mRNA. 
(PDF) 

Figure S3 Metabolic parameters in Bmall ' mice. Body 
weight (A) and daily food intake (B) were measured in 
Bmall — I— and control mice at 14 week of age. (C) Bmall ' 
(red fine) and control mice (black hne) were fed only during the 
night for one week before the sacrifice. Plasma glucose levels were 
measured at the indicated time points (n = 3-6). 
(PDF) 

Figure S4 Phases distribution of SREBPl target genes in 
Bmall~'~ mice. The graph shows the smoothing of phase 
distributions of the genes belonging to the three clusters (green 
line for Al, blue line for A2, magenta line for A3). Only genes 
showing a P-value<0.05 in wild-type mice are plotted. 
(PDF) 

Table SI Sequencing data for the ChlP-seq of SREBPl: 
Number of sequenced and non-redundant tags at each time point. 
(PDF) 

Table S2 Enrichment of motifs for SREBPl, SPl, NFY and 
HNF4 in SREBPl peaks. To estimate the empirical P- Value for 
motifs discovered by MEME, we randomly selected 1000 groups 
of 236 regions and counted the number of matches of the 
indicated motifs in these regions. The random regions were 
selected to bear the same characteristics as the cluster A regions: 
their size and distance to TSS are similar to the ones of cluster A 
regions, and they are close to the TSS of genes expressed in our 
dataset (expression level above the median expression level of eiU 
the transcripts). For all 1000 random groups, we found less 
matches than in the regions associated to SREBP, showing that 
these motifs are enriched with an empirical P-Value<0.001. 
(PDF) 

Table S3 Functional annotation clustering of putative SREBPl 
targets using DAVID tools. The annotation links the sites 
belonging to cluster A to the closest gene irrespective of the 
distance. In total, 219 out of 236 sites have functional annotation. 
The genes associated to each functional term that was significandy 
enriched are shown. P-value = Modified Fisher's Exact Test, 
B = Benjamini corrected P-value. 
(TXT) 
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Table S4 Functional annotation clustering of putati\c SREBP 1 
targets displaying a different combination of binding sites for SPl, 
NFY and/ or HNF4. Functional annotation clustering for the four 
sub-groups of putative SREBP 1 target genes presenting a different 
combination of binding sites for NFY, SPl and/or HNF4. The 
genes associated to each functional term that was significandy 
enriched are shown. P-value = Modified Fisher's Exact Test, 
B = Benjamini corrected P-value. 
PCLSX) 

Table S5 Gene expression profiles in liver tissue of wild type and 
Bmall ' mice for genes reported in Figure 3A and Figure 5 A. 
Microarray gene expression analyses were performed in liver 
samples from wild type (CCNRC samples) and Bmall ' mice 
(NRC samples). Mice were fed only during the night (ZT12— 
ZT24) for one week before collecting liver. Reported data are from 
2 pools/condition (n = 3-5). 
(TXT) 

Table S6 Frequency of HNF4 motifs in putative SREBP 1 
targets with a different temporal expression profile. The 
percentages refer to the frequency of the HNF4 motif that was 
discovered by MEME in the regions under SREBP 1 peaks of the 

three cluster of genes with different temporal expression profile. 
HNF4 motifs were significantly enriched in the peaks of genes 
expressed between ZT4.4 and ZT13 (P-value<0.02 according to 
Fisher's exact test). 
(PDF) 

Table S7 Primer sequences used in qPCR analysis of gene 

expression. 

(PDF) 

Table S8 Primer sequences used in ChlP-qPCR experiments. 
(PDF) 
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